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In the triple ionization of the Li ground state by single photon absorption the three electrons 
escape to the continuum mainly through two collision sequences with individual collisions separated 
by time intervals on the attosecond scale. We investigate the traces of these two collision sequences 
in the classical probability densities. We show that each collision sequence has characteristic phase 
space properties which distinguish it from the other. Classical probability densities are the closest 
analog to quantum mechanical densities allowing our results to be directly compared to quantum 
mechanical results. 
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I. INTRODUCTION 

The theoretical treatment of multiple ionization pro- 
cesses by single photon absorption is highly complex with 
no analytic solution. In the energy domain, the difficulty 
is that one has to account for the correlated motion of the 
electrons in the asymptotic form of the final continuum 
state. In the time domain, this difficulty can be avoided 
at the expense of propagating the fully coupled few-body 
Coulomb problem in time. 

One can surmount the obstacles in the theoretical 
treatment of the triple photo-ionization from the ground 
state of Lithium, for a wide range of energies, by formu- 
lating the four-body break-up process quasiclassically [l| . 
This implies classical propagation of the Coulomb four- 
body problem using the classical trajectory Monte Carlo 
(CTMC) phase space method. CTMC has often been 
used to describe break-up processes induced by particle 
impact 0, EL IE[ with implementations differing usually 
in the way the phase space distribution of the initial state 
is constructed. We use a Wigner transform of the initial 
quantum wave function for the initial state, and this is 
why we call our approach "quasi" -classical. Naturally, 
the electron-electron interaction is treated to all orders 
in the propagation, and any difficulties with electron cor- 
relation in the final state are absent, since the method is 
explicitly time-dependent. The results from the quasi- 
classical formulation for a wide range of energies |1] were 
found to be in very good agreement with experimental re- 
sults [(| 0] as well as theoretical ones available for higher 
excess energies [!, [i[ . 

Moreover, our classical results allow for a detailed anal- 
ysis of the physical processes in terms of the classical 
trajectories: As we have demonstrated, the triply photo- 
ionizing trajectories can be organized in groups according 
to the respective sequence of electron-electron collisions 
According to this collision scheme we have iden- 
tified two main sequences that lead to triple ionization 
from the Li ground state. An indirect verification of the 
collision scheme could be achieved by measurement of 



the electronic angular correlation probability: we have 
shown for excess energies close to threshold [10], that the 
classification scheme of ionizing trajectories can explain 
the electronic angular correlation probability in terms of 
the dominant "T-shaped" pattern of the three escaping 
electrons. The electronic angular correlation probability 
is not yet known experimentally. However, it should be 
measurable with state of the art experimental techniques. 

In the current paper, we explore the manifestations in 
classical probability densities of the two main collision 
sequences the three electrons follow to ionize from the 
ground state of Li. While our previous treatment of the 
collision sequences was on the level of single trajectories 
p], Qij we now treat them on the level of ensemble av- 
erages. Our motivation for doing so is that probability 
densities are the closest classical analog to quantum me- 
chanics. In quantum mechanics the probability density is 
defined directly through the quantum mechanical wave- 
function. In classical physics probability densities can be 
easily computed allowing for a visualization of the dif- 
ferences between classical and quantum mechanical ob- 
servables and for exploring the correspondence principle. 
Our ideas should be a useful tool for identifying and un- 
derstanding collision mechanisms in other systems where 
collision processes play an important role, e.g., strongly 
driven atomic systems. 

Finally, the collision sequences in triple ionization of 
Li take place on an attosecond time scale. While the 
first collision, in each ionization path, occurs around a 
couple of attoseconds after photoabsorption the second 
collision takes place around 70 attoseconds. This is an- 
other interesting aspect of our work: the emergence of 
attosecond laser pulses represents one of the most excit- 
ing develo pm ents in ultrafast laser science over the past 
few years [ll|, LL2| . Attosecond pulses raise the prospect 
of studying electronic wave-packet motion on the time 
scales at which this motion occurs in nature, namely, the 
atomic unit of time (1 a.u.=24 attoseconds). These time 
scales show why attosecond pulses are new tools for ex- 
ploring electronic processes at their natural time-scale 
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and at dimensions shorter than even atomic-dimensions. 
The advancement of ultra-short laser science and its so 
far success in exploring and controlling electronic motion 
fl3^ renders a future direct experimental verification of 
our collision sequences possible. 



II. TIME-DEPENDENT QUASICLASSICAL 
DESCRIPTION OF IONIZATION 

Typically, there are two structurally different con- 
tributions in quantum mechanical matrix elements 
(ipf\0(t)\ijji), the wavefunctions ipij and the operator 
0(t). We call our description quasiclassical because we 
use the full wavefunctions ipi - exactly translated to a 
phase space density p(pf) through a Wigner transform - 
while the subsequent propagation of the density in time 
is performed fully classically. Here, we are interested in 
final states with all three electrons of Li in the contin- 
uum. Since we propagate the entire four body system 
over very long times, we can simply project onto mo- 
mentum states (corresponding to a measurement at the 
detector). In practise, this is done by binning final mo- 
menta of trajectories, very similarly as in the experiment. 



A. The initial phase space density for single 
photon multiple ionization 

The construction of the initial phase space density p(pf) 
in our quasiclassical formulation of the triple photoion- 
ization of Li has been detailed in [l[ , here we give only a 
brief summary. We formulate the triple photoionization 
process from the Li ground state (Is 2 2s) as a two step 
process [H, [3, EH • First, one electron absorbs the photon 
(photo-electron) at time t = t a b s = 0. Then, due to the 
electronic correlations, redistribution of the energy takes 
place resulting in three electrons escaping to the contin- 
uum. It is the latter step that we describe in our formu- 
lation. We first assume that the photo-electron is a ls- 
electron. It absorbs the photon at the nucleus (1*1 = 0), 
an approximation that becomes exact in the limit of high 
photon energy [l6[. The photon could also be absorbed 
by the Li 2s-electron. However, the cross section for pho- 
ton absorption from a Is orbital is much larger than from 
a 2s orbital [l7| • Hence, we can safely assume that the 
photo-electron is a Is electron which significantly reduces 
the initial phase space to be sampled. Also, by virtue of 
their different character the electrons become practically 
distinguishable and allow us to neglect antisymmetriza- 
tion of the initial state. We denote the photo-electron by 
1, the other Is electron by 2 and the 2s electron by 3. Im- 
mediately after photon absorption, we model the initial 
phase space distribution of the remaining two electrons, 
Is and 2s, by the Wigner transform of the corresponding 
initial wavefunction ^(ri = 0,r2,r3), where are the 
electron vectors starting at the nucleus. We approximate 
the initial wavefunction as a simple product of hydro- 



genic orbitals ^(r^) with effective charges Z^ to facil- 
itate the Wigner transformation. The Z{ are chosen to 
reproduce the known ionization potentials i^, namely for 
the 2s electron Z 3 = 1.259 (I 3 = 0.198 a.u.) and for the 
Is electron Z 2 = 2.358 (I 2 = 2.780 a.u.). (We use atomic 
units throughout the paper if not stated otherwise.) The 
excess energy, E, is given by E = — I with E u the 
photon energy and / = 7.478 a.u. the Li triple ionization 
threshold energy. Following these considerations, the ini- 
tial phase space density is given by 

p( 1 ) = yKS(r 1 )d(e 1 + h-u J ) ]J (n, Pi )5( £l + h) 

i=2,3 

(1) 

with normalization constant JV . 

To determine which fraction of ^(7) leads to triple ion- 
ization, the phase space distribution must be followed in 
time. 



B. The evolution of classical phase space densities 

The evolution of a classical phase space density is de- 
termined by the classical Liouville equation, which may 
be obtained within the quantum mechanical phase space 
Wigner formalism [18] by taking the limit h = pH l20j], 

(2) 

The initial phase space values are 

r(o) = 7 , (3) 

and J? c i is the classical Liouville operator which is defined 
by the Poisson bracket {H, }, with H the Hamiltonian of 
the system. In our case H is the full Coulomb four-body 
Hamiltonian. In practice, Eq. ([2|) amounts to discretizing 
the initial phase space, assigning weights to each discrete 
point 7j = (pj (0), qj(0)) according to p(7j), and evolv- 
ing in time each initial condition jj with the Coulomb 
four-body Hamiltonian. This amounts to propagating 
electron trajectories using the classical equations of mo- 
tion (CTMC). Regularized coordinates [2l| are used to 
avoid problems with electron trajectories starting at the 
nucleus. 



C. Different triple ionizing collision sequences and 
their phase space ensembles 

An important finding of our previous studies [TO] is 
that the triple ionizing trajectories can be organized in 
classes according to their ionization-driven properties. In 
particular, we found two main classes consisting of those 
trajectories that triply ionize through the (12,13) colli- 
sion sequence and those that ionize through the (12,23) 
collision sequence. In the first class the path to ioniza- 
tion proceeds with photo electron 1 knocking out, succes- 
sively, electrons 2 and 3. In the second class the photo- 
electron 1 first knocks out electron 2 and then, electron 
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2 knocks out electron 3. More abstractly speaking, each 
class defines an ensemble of trajectories which we label 
a = I and a = II for the (12,13) and the (12,23) collision 
sequences, respectively. 

For completeness, we briefly describe what we define 
as a momentum transferring electron-electron collision 
along a trajectory with time dependent electron posi- 
tions ri(t), i = 1,2,3 (see [lQ|). The term responsible 
for momentum transfer between electrons i and j is their 
Coulomb repulsion V(r^) = r^ 1 , = r^ — rj. Hence, we 
identify a collision between electron i and j (ij) through 
the momentum transfer 

T>ij := - / 2 W(ry) dt = I 2 Fij dt (4) 

under the condition that V(r^(tfc)), fc = 1,2 are local 
minima in time with t2 > ti, while r^- = | r ^ — r j I - This 
automatically ensures that the integral of Eq. (j4j) includes 
the "collision" with a local maximum of V(r^(t)) at a 
time t\ < tM < t2- During the time interval t\ < t < £2, 
all four particles interact with each other. Hence, the 
definition Eq. (|4|) is only meaningful if the collision redis- 
tributes energy dominantly within the subsystem given 
by the three-body Li + -Hamiltonian, H^, of the nucleus 
and the electrons i and j involved in the actual collision, 

Hij =Hi + Hj + 1 , (5) 

with 

dH- ■ 

Hij = » for h < t < t 2 , (6) 

where 

Hi = pi /2 -Z In (7) 

are hydrogenic two-body Hamiltonians with charge Z = 

3 of the Lithium ion. 

Our goal is to investigate whether the two main ion- 
ization sequences we have previously identified using 
Eq. manifest themselves on an ensemble average 
level with properties that clearly distinguish one from 
the other, thus reinforcing the validity of our classifica- 
tion scheme. To this end we need classical observables 
defined over arbitrary phase space ensembles a, in our 
case the two ensembles I and II. 



D. Classical probability densities for observables 
over a classical phase space ensemble 

The probability density ^ a (a, t) to find the value a for 
the observable A at time t under the ensemble a is given 

by 

& a {a, t)= f 6(a - A( 7 , i))p( 7 )d 7 , (8) 



where j a denotes integration over initial phase space 
which contains only those trajectories that belong to the 
ensemble a. The propagation begins at the time £ a b s = 
of photoabsorption. Eq. (|5J) amounts to a) propagating 
all the trajectories of the ensemble a from time £ a b s up to 
t, b) computing for each trajectory the observable A(t), 
c) selecting only those trajectories which satisfy A(t) = a 
and adding together their weights. Note, that the proba- 
bility to find at time t the value a for the observable A(t) 
is given by ^ a (a,t)da. 

Finally, the classical avera ge o f the observable A(t) 
over the ensemble a is simply [22[ 

(A(t)) a = j A( 7 ,^( 7 )d 7 . (9) 

J a 

III. ENSEMBLE AVERAGES OF ENERGY AND 
THE MAIN COLLISION SEQUENCES 

That the ensembles I and II, defined by the two main 
collision sequences, leave different traces on classical av- 
erages is obvious from Fig. [TJ This and all other results 
have been obtained at an excess energy of E = 0.9 eV, 
that is fairly close to threshold [l|. Firstly, we see in 
Fig. that the electron pair potential energies (l/r^j 
(solid line) and {l/r\z) l (dashed line) have well defined 
maxima at t\2 — 1.6 and = 65 attoseconds, even 
when averaged over all trajectories that ionize through 
the (12,13) collision sequence. Note that (1/7*23)! (dot- 
ted line) decreases monotonically. Thus, all trajectories 
triply ionizing through the (12, 13) collision sequence, 
satisfy as an ensemble the first criterion of this sequence, 
namely, maxima in the potential energies of the electron 
pairs 12 and 13 participating in collisions. They also sat- 
isfy as an ensemble the second criterion, that is, while 
the (1/7*12)! potential energy changes during the 12 col- 
lision the energy of the three-body Hamiltonian (Hi2) v 
see Eqs. [5j [6j remains constant and while the (1/7*13) j 
potential energy changes during the 13 collision the en- 
ergy of the three-body Hamiltonian (His) l remains con- 
stant. This is clearly demonstrated in Fig. [TJd: (#12)1 
(thin solid line) remains almost constant near t±2 = 1.6 
as, while it changes around £13 = 65 as, as should be the 
case since during the 13 collision it is the (His) l energy 
(thin dashed line) that is conserved. Similarly, (Hi3)i re " 
mains constant around t\s = 65 as, while it changes near 
ti2 = 1.6 as, since during the 12 collision it is the (#12)1 
energy that is conserved. We also plot (#23)1 (thin dot- 
ted line) to demonstrate that it changes both, around the 
ti2 and t\s collision times. One may summarize the two 
criteria for triple ionizing ensembles of trajectories as 

(A) A maximum in time of the ensemble average 
(l/rij) a defines a collision and its time Uj between 
electrons i and j. 

(B) Near the time of collision t^ the corresponding 
three-body energy (Hij) a of the ensemble remains 
approximately constant. 
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FIG. 1: Averages over ensembles a — I (left, a)-b)) and a — 
II (right, c)-d)), see Eq. ©• The upper panels show the 
interelectronic repulsions, (l/m)^ (solid), (1/7*13) a (dashed), 
and (l/r23) a (dotted). The lower panels show in addition 
the three-body energies with thin lines {H\2) a (solid), (H13) a 
(dashed), {H 2 z) a (dotted). 



These criteria apply also to ensemble II for the collision 
sequence (12,23) as one can see in Figs. [It and[T]i. The 
collisions 12 and 23 are well defined and take place with 
maxima in (1/7*12)11 at t\2 = 1-7 as and in (1/7*23)11 
t 2 s = 69 as, the corresponding three-body energies re- 
main almost constant around these times. 

It is important to keep in mind that the collisions de- 
scribed above are not binary but three-body collisions in- 
volving two electrons and the nucleus. If we were dealing 
with binary collisions then {pi/ 2 + Pj/2 + V | r * — r j|) 
instead of (Hij) would be constant during the ij colli- 
sion, which is not the case as one can easily show. For 
simplicity we will identify in the following a collision by 
ij. However, it is always understood that the nucleus is 
part of the collision and that ij refers to the three-body 
Hamiltonian H^ as given in Eq. ([5]). 



IV. PROBABILITY DENSITIES 

A. Position 

The probability densities for the Cartesian positions 
of each of the three electrons z, given by Eq. (|5J) where 
A(j : t) = qi(jit) with q = x,y or z, do not offer much 
information regarding the identity of the electron pairs 
that participate in each collision, see Fig. [2j However, 
they contain information about the electron pairs dur- 
ing the collisions: Fig. [2] shows that the probability den- 



FIG. 2: Top panel: £P a (xi,t) for electrons i = 1,2, 3 (from 
left to right) for the a — I ensemble; Bottom panel: as for the 
top panel but for £? a (yi,t). 



sity of the Cartesian coordinates of electrons 2 and 3 
remain almost constant until 1.6 and 65 as, respectively, 
the times that electron 1 knocks them out in the 12 and 
13 collisions of the (12, 13) collision sequence. In addi- 
tion, the probability density of the y and z component 
of the photoelectron 1 does not change up to 1.6 as, the 
time of collision of electron 1 with 2. This is expected for 
small times since the photoelectron's initial momentum 
is along the x-axis. Similar conclusions can be drawn for 
the (12,23) collision sequence. 

Due to our choice of initial conditions for the photo- 
electron (electron 1 always starts at r = with initial 
momentum along the x-axis) our model has cylindrical 
symmetry around the x-axis. As a result the probabil- 
ity density of yi is equal to that of Z{. Moreover, since 
the invariance under the parity operation y — > — y of 
the Hamiltonian is not broken through the initial condi- 
tions, Eq. (pQ), the corresponding distributions &(yut) 
are symmetric about yi = 0, the same holds of course 
for the ^-coordinates and for all times, see Fig. [2j These 
symmetry properties extend to observables which respect 
them, such as, e.g., individual electron momenta. 



B. Momentum 

The probability densities of the Cartesian momentum 
components for each of the electrons allow for a deeper in- 
sight into the mechanism of the electron collisions taking 
place during the two main sequences. During an ij col- 
lision the transfer of energy between electrons i and j is 
mediated through their mutual repulsion, V(rij). As dis- 
cussed above, we have defined the collision time tij as the 
time V(rij) = with Vij maximal, see Eq. (j4]). However, 
a collision may formally be defined to last between two 
minima of Vij at times t\ < Uj and £2 > Uj- During the 
time the two electrons approach each other, t\ < t < t^, 



5 



V(rij) > 0, while for Uj <t<t 2 , V(rij) < 0. 

This has different consequences for the correspond- 
ing time evolutions of the individual hydrogenic energies 
Hi(t) of the two electrons, Eq. (|T|), the one which suf- 
fers a net loss of energy in the collision (the impacting 
electron) and the one which overall gains energy in the 
collision (the impacted electron), see Fig. El The latter 
gains energy throughout the collision that is, Hj > for 
t\ < t < t2- On the other hand, the impacting elec- 
tron looses energy up to a time t s , that is Hi < for 
t\ < t < t s with t s > tij (t s > tij follows from Eqs. 
and V(rij) = at time Uj) while for t s < t < the im- 
pacting electron gains energy Hi > 0. From Eqs. ([5]) and 
(|6]) we can determine the change of a hydrogenic energy 
Hj in time during a collision with electron i. On the one 
hand we have 



dlLj 
dt 



Pj ■ Pj 



On the other hand we have 



Pi 



dH 



Inserting Eq. dTTJ) into Eq. (JTDI) leads to 



dt 



= Pj ■ F i. 



(10) 



(11) 



(12) 



which shows, that the change of the hydrogenic energy 
of electron j does not only depend on the modulus of the 
electron-electron force but also on its direction rela- 
tive to the momentum pj of electron j. We recall at this 
point that the present analysis of the collision sequences 
in terms of the rate of change of the single electron ener- 
gies Hi does by no means imply that we have calculated 
the evolution of trajectories with Hi. All observables 
are evaluated with our numerical results for the triple 
ionizing trajectories from the propagated full four-body 
Coulomb Hamiltonian, as we have already pointed out in 
section IIA. 

In Fig. [3] the ensemble averages (Hi) a clearly illustrate 
the difference between the impacting and the impacted 
electron for each three-body collision. For trajectories 
from the ensemble a = I (Fig. [3^i) electron 1 transfers 
energy to electron 2 during the 12 collision as can be seen 
from the sharp decrease of (i^i)i followed by an increase 
beginning at t = 2.4 as while at the same time the energy 
(#2)1 °f the impacted electron 2 increases. The pattern 
is repeated during the 13 collision where a decrease in 
(Hi) 1 is followed by an increase at 79 as, while at the 
same time (H^)i of the impacted electron 3 is increasing. 
The pattern of the hydrogenic energies during collisions 
is also fullfilled for the ensemble II as can be seen in 
Fig. O). 

Describing the collisions using the rate of change of the 
single electron Hamiltonians has the advantage that the 
effect of the nucleus is "folded in" . As a result both, the 
early collision which takes place close to the nucleus and 




FIG. 3: Single electron energy averages (Hi) a of electrons 1 
(solid), 2 (dashed), 3 (dotted) for the ensemble a = I ( a)) 
and a = 11(b)). 



the latter one which takes place far away from it (for both 
ensembles I and II), exhibit exactly the same pattern, see 
Fig. [3l This becomes even clearer when one compares 
Fig. [3] with the probability densities of the momentum 
component along the x-direction for all three electrons in 
Figs. HJ [5] and [6j The momentum along the x-direction 
of the impacting electron 1 in ensemble I switches form 
decreasing to increasing at times 3.5 as for the 12 collision 
and 79 as for the 13 collision. This is a consequence of 
Hi switching sign at 2.4 as for the 12 and at 79 as for 
the 13 collision. The time of 3.5 as, where p Xy i starts 
to increase in the 12 collision, is different from the time 
2.4 as, where Hi switches sign, while both times are the 
same in the 13 collision. The reason is that for the 12 
collision the x, y, ^-coordinates are not equivalent, with 
the transfer of momentum taking place mainly along the 
x-direction, while for the later 13 collision the x, z- 
coordinates are almost equivalent. This is illustrated in 
Fig. [3 with pi • F 2 i = at 2.4 as and p X: i • F 2 i = 
at 3.5 as. The nucleus has a significant effect on the 

12 collision while it has a small one on the 13 collision 
as seen by the more prominent increase of p Xj i at 79 as 
when compared to its increase at 3.5 as. The change 
with time of p Xj i during the 12 collision is due to F21 
and the —dHi/dri force from the nucleus, while in the 

13 collision the change of p Xy i is mainly due to F31. In 
Figs. [5] and [6] we see that p Xj 2 and p x ^ increase during 
the time the respective electrons 2 and 3 are impacted 
by electron 1, in agreement with H2 > and H3 > 
during the 12 and 13 collisions. Similar conclusions can 
be drawn for the a = 11 ensemble. 



C. Inter-electronic angles 

Finally, we discuss the time evolution of the inter- 
electronic angles. For large times and quasi- free motion 

ex p^, the inter-electronic angles refer to the relation 
between positions as well as momenta of the electrons. 
The dynamics in the angle is governed by two principles: 

(A) Collisions between two electrons lead to a minimum 
of the angle Oij between the participating electrons 
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FIG. 4: Momentum distributions & a (p x ,t) for electron 1 for 
the a — I (top) and a = II (bottom) ensemble. The left 
panels show the evolution for short times in greater detail. 
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FIG. 5: Same as Fig. |4] but for electron 2. 
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FIG. 6: Momentum distribution ^ a (p x ,t) for electron 3 for 
the a = I (left) and a = II (right) ensemble. 
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FIG. 7: Averaged angles as a function of time over the a = 
I ensemble. Solid: angle <j> between F21 and pi in a) and 
between F31 and pi in b); Dashed: angle <j) between F21 and 
Paj,i in a) and between F31 and p x ,i in b); Dash-dotted: angle 
4> between F12 and P2 in a) and between F13 and P3 in b); 
Dotted: angle <j> between F12 and p x ,2 in a) and between F13 
and pec, 3 in b). 

i and j, i.e., Oijfaj) « 0, if the collision happens at 
time tij. 

(B) Electrons tend to move away from each other min- 
imizing their mutual repulsive interaction. This 
leads to an interelectronic angle of 180°, if none 
of the electrons suffers a collision through the third 
electron. 

With these two principles, we recognize in Fig. [8] the first 
collision early on (small angle, criterion A) and we also 
infer that a second collision happens around 65 as, but 
not between electrons 1 and 2, since suddenly the in- 
crease of their mutual angle towards 180° (criterion B) is 
stopped and #12 shrinks again towards its final value of 
90°, giving rise to the "T-shape" structure of the three 
escaping electrons [lOj . This is true for both collision se- 
quences, (12,13) and (12,23). The first collision happens 
in both cases between electrons 1 and 2, and then elec- 
tron 3 imposes a second collision with one of the partners 
forming the angle in Fig. [51 namely with electron 1 (up- 
per panels) and electron 2 (lower panels). Finally, since 
in both cases the last colliding electron pair is not the 12, 
612 approaches 90°. 

The evolution of «0*(0i 3 ), Fig. H and ^(<9 23 ), Fig. QUI 
differs much more for the respective two sequences (upper 
and lower panels). However, there is a similarity across 
the two figures, namely the pattern in the upper (lower) 
panel of Fig.[9]is similar to that in the lower(upper) panel 
of Fig. M 

The reason is that in the case of #13 only the a = I en- 
semble (upper panel of Fig. [9J leaves the visible imprint 
of a collision, bringing #13 close to zero while it rapidly 
approaches 180° afterwards, since it is the last collision 
and electrons 1 and 3 move away from each other after- 
wards. The same, but now for electrons 2 and 3, is true 
for the a = II ensemble (lower panel of Fig. fTQj) . 

In the other two panels (lower panel of Fig. [9] and upper 
panel of Fig. [TDJ) one recognizes with the sudden turn 
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away from 180° for the respective interelectronic angle a 
collision with the third electron (criterion B), in case of 
Fig. [9] the collision partner is electron 2 and in case of 
Fig. [TOl electron 1. 

Finally, it is worthwhile to note and understand the 
great variation in the width of the initial distribution for 
the angles across the three figures. ^(#12, 0) is most 
strongly confined to values around zero since the 12- 
collision happens at an early time and electron 1, having 
absorbed the photon energy, starts with relatively high 
velocity close to the origin (position of the nucleus). In 
the short time (1.7 as) before the collision the momen- 
tum vector of electron 1 cannot change substantially, so 
electron 1 keeps its direction. 

The condition for the first collision, r 1(^12) « ^2(^12), 
implies also that the angles #13 and #23 should be similar 
at early times. This is indeed the case, comparing the up- 
per left panels of Figs. I9l and [TOl where in both cases the 
maximum of the initial distribution is around 50°, while 
for the lower left panels the widely spread initial distribu- 
tion is centered about 90°. The latter indicates no clear 
preference in the initial mutual angle between electrons 2 
and 3 reflecting the expectation value of an uncorrelated 
(product) wavefunction for the Lithium ground state as 
used here. Moreover, one should keep in mind that when- 
ever electron 3 is involved one would expect a wider dis- 
tribution due to the larger size of the 2s initial electron 
density compared to the Is density for electrons 1 and 2. 




FIG. 8: Same as Fig. [4] but for the probability density of the 
inter-electronic angle #12, &a(0i2,t). The arrows indicate the 
time of the collision tn. 



CONCLUSIONS 



We have investigated the two main collision paths the 
three electrons follow to escape to the continuum from 
the ground state of Li after single photon absorption on 
an ensemble level. Studying the classical probability den- 
sities for the two ensembles of trajectories correspond- 
ing to the two main collision sequences we were able to 
identify the traces these sequences leave on the classical 
probability densities. Furthermore, we could show that 
each of the two ensembles has unique manifestations on 
the ensemble average level which clearly distinguish one 
from the other. Being able to distinguish the two main 
attosecond time scale collision sequences on an ensem- 
ble level holds promise for a future direct observation of 
these collision sequences with the advancement of ultra- 
short laser technology. 
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FIG. 9: Same as in Fig. [8] but for the inter-electronic angle 
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